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Abstract. In this paper we give, for the first time, a complete description of the late- 
time evolution of non-tilted spatially homogeneous cosmologies of Bianchi type VIII. 
The source is assumed to be a perfect fluid with equation of state p = (7 — where 
7 is a constant which satisfies 1 ^ 7 ^ 2. Using the orthonormal frame formalism 
and Hubble-normalized variables, we rigorously establish the limiting behaviour of the 
models at late times, and give asymptotic expansions for the key physical variables. 

The main result is that asymptotic self-similarity breaking occurs, and is accom- 
panied by the phenomenon of Weyl curvature dominance, characterized by the 
divergence of the Hubble-normalized Weyl curvature at late times. 



PACS numbers: 0420H, 0440N, 9880H 

1. Introduction 

A long term goal in theoretical cosmology is to understand the structure and properties 
of the space of all cosmological solutions of the Einstein field equations (EFEs), with a 
view to shedding light on the evolution of the physical universe from the Planck time 
onward. In particular one wants to study deviations from the familiar Friedmann- 
Lemaitre (FL) models, which describe a universe that is exactly homogeneous and 
isotropic on a suitably large scale. In working towards this goal one makes use of a 
symmetry-based hierarchy of cosmological models of increasing complexity, starting 
with the familiar FL models: 

i) FL cosmologies 

ii) non-tilted spatially homogeneous (SH) cosmologies 

iii) tilted SH cosmologies 

iv) G2 cosmologies 

v) Gi cosmologies 

vi) generic cosmologies 

The terminology used in this hierarchy has the following meaning. A SH cosmology 
is said to be tilted if the fluid velocity vector is not orthogonal to the group orbits, 
otherwise the model is said to be non-tilted. A G2 cosmology admits a local two- 
parameter Abelian group of isometrics with spacelike orbits, permitting one degree of 
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freedom as regards spatial inhomogeneity, while a Gi cosmology admits one spacelike 
Killing vector field. 

An important mathematical link between the various classes in the hierarchy is 
provided by the idea of representing the physical state of a cosmological model at an 
instant of time by a point in a state space, which is finite dimensional for classes i)-iii) 
and infinite dimensional otherwise. The EFEs are formulated as first order evolution 
equations, and the evolution of a cosmological model is represented by an orbit (i.e. 
a solution curve) of the evolution equations in the state space. The state space of a 
particular class in the hierarchy is a subset of the state spaces of the more general 
classes, which implies that the particular models are represented as special cases of 
the more general models. This structure opens the possibility that the evolution of 
a model in one class may be approximated, over some time interval, by a model in a 
more special class. 

Detailed information about the evolution of cosmological models more general 
than FL models can only be obtained using numerical simulations or perturbation 
theory. On the other hand, by introducing suitably normalized variables, one can 

hope to use methods from the theory of dynamical systems to obtain qualitative 
information about the asymptotic regimes of cosmological models, namely 

a) the approach to the initial singularity, characterized by f ^ 0+ , and 

b) the late-time evolution, characterized by ^ — > +00, 

where £ is the overall length scale. Prom a dynamical systems point of view, the 

evolution in the asymptotic regimes is governed by the dynamics on the past attractor 
and the future attractor, respectively. From a physical point of view the asymptotic 
regime ^ — > 0+ corresponds to the approach to the Planck time, while the asymptotic 
regime £ — > +00 could describe the later stages of a particular epoch in the evolution 
of the universe, for example the radiation-dominated epoch in the early universe. 

The above comments provide the background for the present paper, which deals 
with one of the unsolved problems concerning non-tilted SH cosmologies. There are 
two generic classes of non-tilted ever-expanding SH cosmologies, namely Bianchi type 
VIII and exceptional Bianchi type VP ^^g, neither of which have been fully analyzed. 
In this paper our interest lies in the late-time asymptotic regime of Bianchi VIII 
models. Two important results concerning this regime have been proved. Firstly, a 
theorem of Wald (1983) shows that Bianchi VIII models with a cosmological constant 
are asymptotic at late times to the de Sitter model, subject to rather weak restrictions 
on the stress-energy tensor. Secondly, Ringstrom (2001) has recently determined 
rigorously the late-time behaviour of vacuum Bianchi VIII modelsf. We will comment 
on his results later. On the other hand, little is known about non- vacuum models with 
zero cosmological constant, apart from a brief heuristic discussion of models containing 
dust or radiation by Doroshkevich et al (1973) and by Lukash (1974). We likewise 
comment on these results later. In the present paper we give a rigorous analysis of the 
asymptotic dynamics at late times of Bianchi VIII cosmologies whose matter content 
is a perfect fluid with equation of state p = (7 — l)/i, where 7 is a constant. This 
equation of state includes the physically important cases of dust (7=1) and radiation 
(7=1)- 

To achieve our goal we formulate the EFEs for the Bianchi VIII cosmologies 
as an asymptotically autonomous system of ordinary differential equations, using 
the Hubble-normalized variables first introduced by Wainwright and Hsu (1989) and 

t Barrow and Caspar (2001) have also studied this problem but their results are not conclusive. 
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modified by Wainwright et al (1999)t and Nilsson et al (2000) in their study of the 
more special Bianchi VIIq models. We are then able to apply a theorem of Strauss 
and Yorke (1967) concerning the solutions of asymptotically autonomous systems of 
differential equations. 

The paper is organized as follows. In section 2 we present the evolution equations 
for SH cosmologies of Bianchi type VIII using the orthonormal frame formalism and 
Hubble- normalized variables. The equations are derived in detail in Wainwright 
and Ellis (1997)| (sec chapters 5 and 6). A change of variable then leads to a 
new form of the evolution equations that is adapted to the oscillatory nature of the 
Bianchi VIII models. In section 3 we present the main results, namely theorem 3.1 
and corollary 3.1, which give the limits in the late-time regime of the expansion- 
normalized variables and of certain physical dimensionless scalars, thereby describing 
the asymptotic dynamics. In section 4 we strengthen the results of the theorem and 
corollary by giving the asymptotic form of the expansion-normalized variables and 
physical scalars at late times. The notions of asymptotic self-similarity breaking 
and Weyl curvature dominance are discussed in section 5. Finally, in section 6 we 
conclude by giving an overview of the asymptotic dynamics of non-tilted SH perfect 
fluid cosmologies, noting that the present paper and the accompanying paper Hewitt 
et al (2002) fill the gaps that exist in the literature. 

There are four appendices. Appendix A contains the proof of the fact that Bianchi 
VIII universes are not asymptotically self-similar at late times. Appendix B contains 
the proof of theorem 3.1 while in appendix C we provide an outline of the derivation 
of the asymptotic forms at late times. Finally, in appendix D we give expressions for 
the components of the Weyl curvature tensor in terms of the expansion-normalized 
variables. 

2. Evolution equations 

In this section we give the evolution equations for non-tilted SH cosmologies of Bianchi 
type VIII. As described in WE (see pp 124-5), we use expansion- normalized variables 



defined relative to a group-invariant orthonormal frame {eg}, with Cq = u, the fluid 
4- velocity, which is normal to the group orbits. 

The variables S± describe the shear of the fluid congruence, and the Na, 
a = 1,2,3 describe the intrinsic curvature of the group orbits. The models of Bianchi 
type VIII are described by the inequality N1N2N3 < 0. Without loss of generality, we 
assume 



(E+,E_,7Vi,7V2,A^3), 



(2.1) 



^1 <0, 



N2 > 0, 



N3 > 0. 



(2.2) 



It is convenient to define 



N+^^(N2 + Ns), 

and replace (2.1) by the state vector 
(I]+,S_,7Vi,7V+,iV_). 



(2.3) 



(2.4) 



t This work will henceforth be referred to as WHU. 
I This work will henceforth be referred to as WE. 
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The restrictions (2.2) become 

Ni < 0, Nl- 3N^ > 0, N+> 0. (2.5) 

The state variables (2.1) and (2.4) are dimensionless, having been normalized 
with the Hubble scalarf H, which is related to the overall length scale i by 

H=j, (2.6) 

where the overdot denotes differentiation with respect to clock time along the fluid 
congruence. The state variables depend on a dimensionless time variable r that is 
related to the length scale £ by 

e = £oe^, (2.7) 

where £o is a constant. The dimensionless time r is related to the clock time t by 

as follows from equations (2.6) and (2.7). In formulating the evolution equations we 
require the deceleration parameter q, defined by 

££ 

.= -p (2.9) 

and the density parameter Q, defined by 

n=^. (2.10) 

In terms of the variables (2.4), the evolution equations (6.9) and (6.10) in WE 
become 

= -(2 - q)^+ ~ 2Nl + iiVi(iVi - iV+), 
S'_ = -(2 - q)Y,_ - N_{2N+ - Ni), 

N[ = (g-4S+)iVi, (2.11) 

= (g + 2S+)Ar+ + 6S_Ar_, 
N'_ = {q + 2S+)iV_ + 2S_iV+, 

where 

q = 2{El + E^) + i(37 - 2)n, (2.12) 

n = l-j:l-^l-Nl- ^Nf - {-^NiN+) . (2.13) 

and ' denotes differentiation with respect to t. For future reference we also note the 
evolution equation for O: 

Q' = [2g - (37 - 2)]0. (2.14) 

The physical requirement > 0, in conjunction with (2.5), implies that the 

variables E+, S_, A'^-, A^i and NiN+ are bounded, but places no restriction on N-i- 
itself. In fact, it will be shown in appendix A (see proposition A.l) that if ^ and 
I ^7^2, then for any initial conditions 

lim N+ = +0O. (2.15) 

r— >+oo 

f On account of (2.6), H is related to the rate of volume expansion of the fluid congruence according 

to H = ie. 
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Inspection of the equations for and N'_ shows that this unboundedness of Nj^ will 

induce rapid oscillations in E_ and A'^_ at late times, which further complicates the 
dynamics. A similar situation occurs in the Bianchi VIIq models analyzed in WHU 
(see p 2583). The first step in analyzing the dynamics at late times (r — > +00) is 
to introduce new variables which arc bounded at late times and which enable us to 
isolate the oscillatory behaviour associated with S_ and A^_. We define 

S_ = i?cosV', N- = Rsmip, 

1 , 1 (2.16) 

N+' 3 +' 

where i? > 0. 

In terms of the new variables (S+, i?, Z, M, tp), the evolution equations (2.11) can 
be shown to have the following form 

= -(2 - g)E+ - i?2 + + SM^Z'^ + (1 + S+)i?2 cos 2^, 

R' = [{Q + S+ - 1) + (ii^ - 1 - S+) cos 2^ - |MZ^ sin2V'] R, 

Z' = [Q - S+ + i?2 (cos 2^ + |M sin 2tl))] Z, ^2 17) 

M' = -[Q + 2S+ + i?2 (cos 2tp + 3M sin 2V;)] M, 



1 



where 



and 



V-' = [2+(l + S+)Msin2V' + |M2z2(l-cos2V')] , 
Q = 2S^ + i?2 + i(37-2)r2, (2.18) 



n = l-i:l- R^ - Z^ - iM'^Z'^. (2.19) 
The evolution equation for f2 becomes 

= [2Q - (87 - 2) + 2i?2 cos 2il)] O. (2.20) 
The restrictions (2.5) arc equivalent to 

ZM^R^ sin2 V < 1, M > 0, i? > 0, Z > Q. (2.21) 

3. Limits at late times 

In this section we present a theorem which gives the limiting behaviour as r ^ +00 of 
non-tilted Bianchi VIII cosmologies when the equation of state parameter 7 satisfies 
1 ^ 7 ^ 2. As a corollary of the theorem, we obtain the limiting behaviour of 
certain dimensionless scalars that describe physical properties of the models, namely 
the density parameter f2, defined by (2.10), the shear parameter S, defined by 

= (3-1) 

where Oab is the rate-of-shear tensor of the fluid congruence, and the Weyl curvature 
parameter W, defined by 

where Eab and Hab are the electric and magnetic parts of the Weyl tensor, respectively 
(see WE, p 19), relative to the fluid congruence. 
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In terms of the Hubble- normalized variables, the shear parameter is given by 

= Y.\+ cos^ V, (3.3) 

which follows from (2.16) in conjunction with equation (6.13) in WE. The formula for 

the Weyl curvature parameter is more complicated and is provided in appendix D. 

The main result concerning the limits of R, Z and M is contained in the 
following theorem. One of the limits depends on requiring that the model is not 
locally rotationally symmctricf (LRS). 

Theorem 3.1. For all non-tilted SH cosmologies of Bianchi type VIII, with equation 
of state parameter 7 subject to 1 ^ 7 ^ 2, the Hubble-normalized state variables 
{T,+ ,R,Z,M) satisfy 

hin (E+, R, Z, M) = (i, 0, ^, 0) . (3.4) 

If the model is not LRS, then 

M 

hin - = 0. (3.5) 

r— >+oo Jrl 

Proof. We first consider models which are not LRS. It follows immediately from 
(2.15) and (2.16) that 

lim M = 0. (3.6) 

T — ^ + 00 

Furthermore, since E+ and Z are bounded, it follows from the if) evolution equation 
in (2.17) that 

lim if} = +00. (3.7) 

r — i-+oo 

The trigonometric functions in the DE (2.17) thus oscillate increasingly rapidly as 
T — > +00. In order to control these oscillations, we define new gravitational variables 
S_|_, R and Z according to 

S+ = E+ - i(l + Y.+)R^M s,m.2i), 

R =i?[l- i(i?2 - 1 -S+)Msin2V'] , (3.8) 

Z = Z [I - \R^ M sai2i}}] , 

motivated by the analysis of WHU. The evolution equations for these "barred" 
variables, which can be derived from (2.17) and (3.8), have the following form 

% = -(2 - g)E+ -R? + Z'' + MBj,^, 

R' ={Q + f:+-l + MB^)R, (3.9) 
Z' ={Q-f:+ + MBz)Z, 

where 

Q = 2S^ + i?2 + 1 (3^ _ 2)(1 - - i?2 „ 2^), (3.10) 

and the B terms are bounded functions in R, Z and in M and tp for r sufficiently 
large. The essential idea is to regard M and ip as arbitrary functions of r subject only 
to (3.6). Thus, (3.9) is a non-autonomous DE for 

x = (S+,^,Z), (3.11) 

t See for example, WE p 22. We note that the LRS Bianchi VIII models are described by the 
invariant subset S_ = N- = 0, equivalently, R = 0. 
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of the form 

x' = f{x)+g{x,T), (3.12) 

where 

g{x, t) = M{t) (Sg^ , RB^, ZBz) , (3.13) 

and the form of f{x) can be read off from the right-hand side of (3.9). Since 
hm g{x,T) = 0, 

r— *+oo 

as follows from (3.6), the DE (3.9) is asymptotically autonomous (see for example, 
Strauss and Yorke 1967). The corresponding autonomous DE is 

X = fix), (3.14) 

where 

x = {t+,KZ). (3.15) 
Using standard methods from the theory of dynamical systems, we first show that 
lim {t+,R,Z)=U,0,^). (3.16) 

r— >+(X) \^ ^ / 

Details are provided in appendix B.l. We then use a theorem from Strauss and Yorke 
(1967) (see theorem B.l in appendix B) to infer that the solutions of (3.12) have the 
same limits as the solutions of (3.14), namely 

Km (S+,i?,Z) = (^i,0, ^) . (3.17) 

Details are provided in appendix B.2. The limit of a; = (S+, i?, Z) follows immediately 
from this result in conjunction with the definitions (3.8). Finally, (3.5) is derived in 
appendix B.3, using (3.4). 

The proof of theorem 3.1 for the case of LRS models is straightforward. Since 
R = 0, the oscillatory terms in the evolution equations (2.17) drop out, and the 
variable ip becomes irrelevant. Since (3.6) still holds, the resulting system of equations 
is a special case of the autonomous DE (3.14), with the result that S+ and Z have 
the same limits as in the non-LRS case. □ 

Comment. Our proof of theorem 3.1 can be adapted to the case of vacuum models, 

in which case we recover the results of Ringstrom (2001) (see theorem 1, p 3793). 

Corollary 3.1. For all non-tilted SH cosmologies of Bianchi type VIII, with equation 
of state parameter 7 subject to 1 ^ 7 ^ 2, the shear parameter S and the density 
parameter Q. satisfy 

lim S = i, lim O = 0. 

T— »+cx) r— *+oo 

The Weyl curvature parameter W satisfies 



and 



lim W = +00, if the model is not LRS, 



lim W = 0, if the model is LRS. 
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Proof. These results follow directly from theorem 3.1 and equations (2.19), (3.3), 
(D.3) and (D.4). In particular, if the model is not LRS, it follows from (D.3) and 
(D.4) that since E_|_, R and Z are bounded, and lim^^+oo M/R = 0, that 



2 R 

yv = — [i + o{M)], 



(3.18) 



as T ^ +00. □ 

We have focussed our attention on models for which the source is a classical fluid, 
i.e. 7 lies in the range 1 < 7 < 2. To conclude this section we statef the limiting 
behaviour of the scalars O, S and W for values of 7 in the range < 7 < 1: 



lim O 



lim S = 

T— *+00 



Um W = •{ 

■— '■+00 




if < 7 < I 



if 



|<7<1 



if 
if I < 7 < 1 



<7<i 



1(37 -2)(1 



7), 



+00, 



if 
if 

if 
if 



0<7< i 

i<7<| 

7=1 
|<7<1 



(3.19) 



(3.20) 



(3.21) 



For the LRS models, we note that (3.19) and (3.20) also hold, while (3.21) becomes 



lim W 



0, 



K37-2)(l-7), 



if 
if 



o<7s; |, 

|<7<1. 



(3.22) 



In the case < 7 < |, these results follow from the fact that for this range of 
7, all Bianchi models are asymptotic at late times to the flat FL model (see WE, 
theorem 8.2, p 174). The case | < 7 < 1 can be treated by the methods described 
in this paper. Note that the models with | < 7 < 1 provide a bridge between the 
inflationary models| with < 7 < | which isotropize, and the models with a classical 
fluid which do not, showing that 7 = | is a bifurcation value. A second bifurcation 
occurs at 7 = I , at which value W makes the transition to unbounded growth. 



4. The asymptotic solution at late times 

In this section we give asymptotic expansions as r ^ +00 for the state variables , 
R, Z, and M that describe the class of Bianchi VIII cosmologies, and for the key 
dimensionless physical scalars fi, E and W. We also give the asymptotic expansion 
for the Hubble scalar H, and the asymptotic relationships between clock time t and 
the dimensionless time r. The angular variable xIj does not play a major role in the 
asymptotic expansions. For brevity, we note its decay rate is governed by 

i^' = ^[1 + 0{M)], 

as T — > +00, as follows from its evolution equation in (2.17). 

f These results were reported earlier in Waiiiwright (2000) (see p 1050). 

X Note that if < 7 < ^, the deceleration parameter q is negative, as follows from equation (2.12). 
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Although the limiting values of the state variables as r — > +00, given in 
thcorcni 3.1, arc valid for all values of the equation of state parameter 7 in the range 
1 ^ 7 ^ 2, it turns out that the asymptotic expansions depend in a significant way 
on 7. In particular, the case 7 = 1 has to be treated separately!. We now give the 
asymptotic expansions in the two cases, 7=1 and 1 < 7 ^ 2, assuming that the 
model is not LRS. In the equations below, Cr, Cm, Co and Ch are positive constants 
andT= is^. 

State variables: 

7=1 1<7<2 

S+ = i - ir-i [1 + 0{T)] S+ = i - ir-i [1 + 0{T)] 

R = Crt-^ [1 + 0{T)] R = ir-i/^ [1 + 0{T)] (4.1) 

^ = - i^^"' [1 + m ^ = f + [1 + 0{T)] 

M = Cmt^I'^Q-^I'^^ [1 + 0{T)\ M = Cmt^I^Q-^I'^^ [1 + C(T)] 



Density parameter: 



\CnT-i/2e-3(7-i)r[i + o(T)], if 1<7<2 



Anisotropy parameters: 



i-ir-Ml + 0(T)], if 7=1 

i- ir-^sin^V [1 + ^5(2^)], if 1<7<2 



^r-5/V/2-[l + 0(T)], 



if 7 = 1 



Hubble scalar: 



Clock time: 



Cj/ri/V3/2"[l + 0(T)], if 1<7<2 



2 ^-i/2e3/2ni + 0(T)], if 7 = 1 



3Ch 
2 



T-i/V/2^[i + 0(r)], if i<7s;2. 



(4.3) 



W = <; (4.4) 

r-''/^e^/^^l + 0{T)], if 1<7<2 



rC^ri/2e-3/2-[l + 0(T)], if 7 = 1 



(4.6) 



Our derivation of the above asymptotic expansions is given in appendix C. While 
not fully rigorous, it is nevertheless quite convincing. We first derive the expansions for 
the solutions x{t) = {T,^,R,Z) of the autonomous DE (3.14), using centre manifold 

t Note that in some parts of the proof of theorem 3.1, the case 7 = 1 has to be treated separately. 
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theory. Secondly, it follows from the proof of theorem 3.1 that M decays to zero 

exponentially. This fact enables us to verify that the expansions are compatible with 
the exact evolution equations (3.9), in the sense that terms cancel at the appropriate 
orders. 

The derivation of the expansion for f7 depends on 7. In the case 7=1, the 
expansion follows immediately from (2.19). The case 1 < 7 ^ 2 is more complicated 
and is treated in appendix C.l. Once Q. is known, the expansion for H can be 
obtained directly from equation (2.10), since it follows from the contracted Bianchi 
identities that the matter density ji is given hy ji = e~'^'^'^ (see for example, WE, 
equation (1.99)). Finally, knowing H, the expansion for the clock time t can be 
obtained from equation (2.8). 



LRS models 

We note that the asymptotic expansions (4.1)-(4.6) were derived subject to the 
assumption that the model is not LRS, i.e. that the variable R is not identically zero. 
A detailed analysis shows that results for LRS models with 7 = 1 can be obtained by 
setting Cn = in equations (4.1)-(4.3) and (4.5)-(4.6). The expansion for the Weyl 
curvature parameter W becomes 

W=\t-^[1 + 0{T)], (4.7) 

For 1 < 7 ^ 2, the asymptotic expansions do not specialize to the LRS case R = 0. 
A detailed analysis shows that the variables approach their limiting values at an 
exponential rate. Specifically, 

S+ = i + 0(e-^-), 

^ ' 4.8 

M =C7Me-3/2^[l + 0(e-^^^)], 

O =C7ne-3(T-i)^[l + 0(e-^^^)], 

where Ai, A2 and A3 are positive constants. It follows that W tends to zero 
exponentially. In conclusion, the main difference is that the LRS models do not 
exhibit Weyl curvature dominance. 



5. Discussion 



The most significant feature of the late-time asymptotic regime of non-tilted SH 
perfect fluid cosmologies of Bianchi type VIII is that the models are not asymptotically 
self-similar, or in other words, the evolution of the models at late times is not 
approximated by a self-similar model. We can draw this conclusion because the orbits 
that describe the models in the Hubble-normalized state space escape to infinity, and 
hence do not approach an equilibrium point (see Wainwright 2000, p 1044; note that 
equilibrium points in the Hubble-normalized state space correspond to self-similar 
Bianchi cosmologies). The Bianchi VIII models share this feature with the Bianchi 
VIIo models (see WHU). However, the Bianchi VIII models differ from the Bianchi 
VIIo models in two important ways. Firstly, the Bianchi VIII models with 1 ^ 7 ^ 2 
are vacuum- dominated at late times, that is, the density parameter ft satisfies 

lim n = 
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(see corollary 3.1). Secondly, the expansion is anisotropic at late times in the sense 
that the shear parameter S satisfies 

lim S = i 

(see corollary 3.1). On the other hand, for the Bianchi VIIq models, the shear 
parameter E tends to zero at late times whenever the equation of state parameter 
satisfies 7 ^ | (see WHU, theorem 2.3). 

The breaking of asymptotic self-similarity in Bianchi VIII and Bianchi VIIq 
cosmologies manifests itself in the behaviour of the Weyl curvature tensor, which 
describes the intrinsic anisotropy in the gravitational field. Both classes of models 
exhibit what has been referred to as Weyl curvature dominance] (sec WHU, p 2588) 
which refers to the fact that Hubble-normalized scalars constructed from the Weyl 
tensor are unbounded. For Bianchi VIII models with 1 ^ 7 ^ 2, Weyl curvature 
dominance manifests itself through the unbounded growth of the Weyl parameter W 
as T ^ +00 (see corollary 3.1), as is also the case for the Bianchi VIIq models with 
1 < 7 < 2 (see WHU, theorem 2.4 and equation (3.40) noting that fi = ^(4-37)). On 
the other hand, for Bianchi VIIq models with 7=1 the Weyl parameter W approaches 
a finite limit and the unboundedness occurs in the Hubble-normalized derivatives of 
the Weyl tensor. 

We summarize the rate of growth of W below, expressed in clock time i, as 
t +00. The results for Bianchi VIIq models follows from equations (3.26), (3.34) 
and (3.40) in WHU. For Bianchi VIIq models, 

L^O, if 7 = 1, 

f2(7-l)/7, if 1<7<|, 



if 7=3' 



(lnt)3/2' 

^(2-7)/7^ if |<7<2, 

while for Bianchi VIII models, 

(^''*^'' ^ (5.2) 

h^' '<^^'- 
We note that rate of growth is largest for Bianchi VIII models with 7 in the range 
1 < 7 ^ 2, and is independent of 7 in this range. For Bianchi VIIq models, the 
maximum rate of growth occurs for the radiation equation of state, but is less than 
that for the Bianchi VIII models. 

We now comment on the relation between the late-time regime of perfect fluid 
Bianchi VIII models and the corresponding vacuum models, as studied by Ringstrom 
(2001). It follows from corollary 3.1 and equation (3.19) that the value 7 = 1 is a 
bifurcation value governing the transition to a vacuum-dominated late-time regime 
for values 7 > 1. This bifurcation manifests itself in the rate at which the density 
parameter O tends to zero. If 7 = 1, decays to zero at a power law rate while if 
1 < 7 < 2, decays exponentially (in terms of r; see equation (4.2)). In other words, 
if 1 < 7 ^ 2, orbits approach the vacuum boundary f7 = in state space exponentially 

f This notion has recently been incorporated into a classification of the late-time dynamics of SH 
models by Barrow and Hervik (2002), who refer to it as extreme Weyl dominance (see section 5). 
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Table 1. References for the asymptotic dynamics of non-tilted SH cosmologies with 
perfect fluidf source and equation of state p = (7 — 



singular regime 



late-time regime 



class A, type I, II, VIq 
class A, type VIIq 
class A, type VIII 
class A, type IX 
non-cxccptional class B 
exceptional class B 



WE, ch 5 
Ringstrom (2000) 
WE, ch 5 
Ringstrom (2000) 
WE, ch 6 

Hewitt et al (2002) 



WE, ch 5 

WHU and Nilsson et al (2000) 
present paper 
not applicable 
WE, ch 6 

Hewitt et al (2002) 



fast, SO that the loading asymptotic behaviour of the perfect fluid models, as given in 
equations (4.1) (4.6), coincides with that of the vacuum models. 

Wc conclude this section by providing a link between our results and the earlier 
analysis of Doroshkevich et al (1973) (see p 742) and Lukash (1974) (see p 168), 
referred to in the introduction. In these papers the authors give an approximate form 
of the line-element at late times for Bianchi VIII models containing dust or radiation 
and of the energy density of the source. In particular, the energy density is 



where t is clock time. Our asymptotic expansions for fi, H and t in equations (4.2), 
(4.5) and (4.6), in conjunction with (2.10), lead to the asymptotic form for /j. which 
agrees with (5.3). 

6. Overview 

In this concluding section we use the results of this paper and of the accompanying 
paper Hewitt et al (2002) to give an overview of the dynamics in the asymptotic 

regimes of non-tiltcd SH perfect fluid cosmologies. These cosmological models can be 
grouped into three main subclasses, following Ellis and MacCallum (1969): 

i) class A models (Bianchi type I, II, VIq, VIIq, VIII and IX), 

ii) non-exceptional class B models (Bianchi type IV, V, VI/j, and Vllh), 

iii) exceptional class B models (Bianchi type Vlh with h = denoted VP ^^g). 

We refer to WE (pp 36-7, 41-2) for a summary of this classification. 

The dimensions of the Hubblc-normalizcd state space for each Bianchi type 
are shown in figure 1. We note that the dimension gives the number of arbitrary 
parameters in the corresponding family of solutions (see Wainwright and Hsu (1989), 
p 1419). 

Table 1 gives the references which contain the most comprehensive descriptions 
of the various classes. All of these papers make use of the so-called Hubble-normalized 

variables within the framework of the orthonormal frame formalism, and apply 
techniques from the theory of dynamical systems. This approach highlights the role 




(5.3) 



t For a complete description of the late-time regime for vacuum class A models, we refer to Ringstrom 
(2001). 




played by self-similar solutions. These works in turn refer to related papers, describe 
other methods, specifically the Hamiltonian approach and the metric approach (see 
WE, pp 5-6), and give information about the historical development. 

The table shows that with the appearance of the present paper and the paper 
Hewitt et al (2002), there is now available a complete description of the dynamics of 
non-tilted SH cosmologies with perfect fluid source, in the two asymptotic regimes. It 
should be noted, however, that some of the results are based on heuristic arguments 
and numerical experiments. We shall mention below what remains to be proved. 

The principal features of the asymptotic dynamics are as follows. Firstly, as 
regards the singular asymptotic regime, the three generic classes, namely Bianchi 
type VIII, IX and Vl^-^^g (see figure 1) have an oscillatory singularity, and the 
asymptotic behaviour is described by a two-dimensional attractor containing the 
Kasner equilibrium points in the Hubble-normalized state space. We note that the 
existence of this attractor has only been proved for the Bianchi IX models (see 
Ringstrom (2000)). For non-generic models the dynamics near the singularity is 
remarkably simple in the sense that the models are asymptotically self-similar, being 
approximated by a Kasner solution. 

Secondly, as regards the late-time regime, the dynamics depends crucially on 
whether or not the Hubblc-normalizcd state space is bounded. For the classes VIIq 
and VIII the state space is unbounded, and the solutions exhibit Weyl curvature 
dominance (see section 5). For the remaining classes, the state space is bounded and 
the models are asymptotically self-similar. In particular, all class B models, including 
the generic class VF ^^g, are asymptotically self-similar, although a complete proof of 
this fact, except in the vacuum subcase, has not been given, due to a lack of success 
in finding a monotone function for the evolution equations. 

t The Bianchi V models, for which the state space is two-dimensional, do not fit in a natural way 
into this diagram. 
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Referring to the hierarchy in the introduction, we conclude by giving some 
suggestions for future research. Firstly, much work remains to be done before 
the asymptotic dynamics of tilted SH models are fully understood. While various 
formulations of the EFEs have been given (see WE, p 175) the only tilted models that 
have been analyzed in detail is the full class of Bianchi II models (see Hewitt et al 
2001) and a special class of Bianchi V models (see Hewitt and Wainwright 1992). The 
difficultly in analyzing the general class of tilted SH models is highlighted by the fact 
that at this time not all of the self-similar solutions have been found. Secondly, it is a 
long-standing conjecture that the occurrence of oscillatory singularities in SH models 
of Bianchi types VIII and IX implies that oscillatory singularities will occur in generic 
spatially inhomogcncous models. Over the past five years, considerable numerical 
and analytical evidence has been provided to support this conjecture (see Weaver et 
al 1998, Berger and Moncrief 1998 and Berger et al 2001). Likewise, the occurence 
of Wcyl curvature dominance in SH models of Bianchi types VIIo and VIII leads us 
to conjecture that this behaviour will also occur in spatially inhomogeneous models, 
perhaps as generic behaviour. The simplest class in which this behaviour could occur is 
the class of G2 cosmologies, since the Bianchi VHq models are contained as a subclass. 
It would thus be of interest to investigate the G2 cosmologies from this point of view. 
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Appendix A. The limit of Ar_|_ at late times 

In this appendix we prove (2.15), concerning the limit of the Hubble- normalized 
variable in the late- time regime. This result is restated as proposition A.l below. 

Proposition A.l. For all non-tilted SH cosmologies of Bianchi type VIII, with 
equation of state parameter 7 subject | ^ 7 ^ 2 and density parameter CI satisfying 
0^0, the Hubble-normalized variable iV_|_ satisfies 

lim iV+ = -\-oo. (A.l) 

T— »-t-00 

Proof. The proof relies on concepts from dynamical systems theory, in particular, 
the notion of w-limit sets (see for example, WE, p 99) and the monotonicity principle 
(WE, theorem 4.12, p 103). 

The Bianchi VIII state space S is defined by the inequalities (2.2) and Cl ^ 0. 
The function 

X={NiN2N:,f, (A.2) 
is positive on S and satisfies 

x' = 6gx, (A.3) 

as follows from (2.3) and (2.11), where q is given by (2.12). If | < 7 < 2 and O > 0, 
q is positive and hence x is increasing along orbits in S. On the other hand, if = 
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or 7 = |, then q may equal zero if S_|_ = S_ = 0. Since S_|_ = S_ = is not an 
invariant set of the cvohition equations (2.11), q cannot remain zero and hence x is 
again increasing along orbits in S. 

In order to apply the monotonicity principle, we consider the set S\S, where S 
is the closure of S. From the definition of S, it follows that S \ S is defined by the 
following restrictions: 

N1N2NS ==0, Ni^ 0, N2 ^0, TVg ^ 0. (A.4) 

By (A. 2), X is defined and equal to zero on S\S. The monotonicity principle implies 
that for any point x G S, the w-limit set iv{x) is contained in the subset of S\S that 
satisfies the condition lim^^s xiv) Oj where s G S\S and y G S. This subset is the 
empty set, since x = on 5 \ 5. Therefore we conclude that lo{x) = </> for all x £ S. 

Finally, suppose that (A.l) does not hold for each orbit in S. Then there exists 
a number 6 > such that for all tq, there exists a t > tq such that N+{t) < b. Since 
the other variables are bounded, then the orbit x{t) lies in a compact set 5* C and 
hence has a limit point in S, contradicting u){x) = (j). Therefore (A.l) holds. □ 



Appendix B. Details of the proof of theorem 3.1 

The proof of theorem 3.1 is based on a result of Strauss and Yorke (1967) (sec corollary 
3.3, p 180) concerning asymptotically autonomous DEs, stated as theorem B.l below. 
Consider a non- autonomous DE 

x' = f{x)+g{x,T), (B.l) 

and the associated autonomous DE 

x' = fix), (B.2) 

where f : D ^ M."', g : D x M. ^ and D is an open subset of M". It is assumed 
that 

Hi : lim g{w{T), r) = for every continuous function w : [tq, +00) D 

r— >+oo 

and 

H2 : any solution of (B.l) with initial condition in D is bounded for t^tq, 
for some tq sufficiently large. 

Theorem B.l. // Hi and H2 are satisfied and any solution of (B.2) with initial 
condition in D satisfies 

lim x{t) = a, 

T— * + 00 

then any solution of (B.l) with initial condition in D satisfies 
lim x{t) = a. 

T — ^ + 00 



We make use of this theorem in appendix B.2. 
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Figure Bl. Orbits in the invariant set S^. 

Appendix B.l. Limits at late times of{T,+ ,R,Z) 
The components of the DE (3.14), x = f{x), are given by 
t'+ = -{2-Q)t+-R^ + Z^, 

it' ={Q + t+-l)R, (B.3) 
Z' ={Q-t+)Z, 

where 

= 2E^ + i?2 + i (37 - 2)f2, (B.4) 

Cl = l-tl- If - Z^. (B.5) 
One can also form an auxihary DE for Ct using (B.3) and (B.5) to find that 

Cl' = [2Q - (37 - 2)]Cl. (B.6) 
We consider the state space S of the DE (B.3) defined by the inequahties 

i? > 0, Z > 0, n>0. (B.7) 

These inequalities in conjunction with (B.5) imply that the state space S is the interior 
of one quarter of a sphere. We also consider the dynamics on the invariant set 
defined by the following restrictions: 

: n = 0, R>0, Z> 0. 

The DE (B.3) admits a positive monotone function 

. = i. (B^8) 
which satisfies 

X' = 3(1 - 7)X, (B.9) 

on the set S. Thus, if 7 ^ 1 there are no equilibrium points, periodic orbits and 

homoclinic orbits in S (sec WE, proposition 4.2). It is immediate upon integrating 
(B.9) and using the boundcdncss of R and Z that for any x G S 

Lo{x) C 5o, if 1 < 7 ^ 2. (B.IO) 
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Figure B2. Orbits in the invariant set = 2k, 7 = 1. 

RZ 



The flow on the invariant set S^-^ is depicted in fignre Bl, which shows the 
projection of the surface = onto the S_|_i?-plane. The essential features are the 
existence of three equilibrium points 

K±: (S+,4Z) = (±1,0,0), 



V : 



E,Z)=(i,0,^) 



which lie on the boundary of S'jj, and the fact that there are no periodic orbits or 
flow-connected heteroclinic cycles on S^. Thus the only potential w-limit sets in 
are the equilibrium points K='= and V and hence for any x G S, the w-limit set is one 
of K± or V. The point K+ can be excluded since it is a local source in S. The point 
K~ can be excluded by considering the evolution equation for Z, which is of the form 

Z' = z{±+,R, Z)Z. 



Since z{\i ) = z{—l, 0, 0) = 3 and Z = at K , it follows that Hm^ 



Z ^ and 



hence that an orbit in S cannot be future asymptotic to K 
u}{x) = V for any x G S. Equivalently, 



We thus conclude that 



lim 

•— >+CX) 



f 1 ^ 



Finally, we consider the case 7=1. 
describes a conserved quantity 



RZ 



(B.ll) 

By (B.9) the function x defined in (B.8) 

(B.12) 



where fc > is a constant that depends on the initial conditions. We see that for 
all fc > the surfaces described by (B.12) foliate the state space S and intersect the 
vacuum boundary $1 = at i? = and Z = {] (sec figure B2). 

We can use similar techniques from before to conclude that the equilibrium points 
K='= are asymptotically unstable. Since the orbits are constrained to lie on the two- 
dimensional family of invariant sets defined by (B.12), it follows that ijj{x) = V for 
any x £ S. Thus any solution of the DE (B.3) in the set S also satisfies (B.ll) if 
7 = 1. 
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Appendix B.2. Limits at late times of (S+, ^, Z) 
We now apply theorem B.l to prove (3.17), namely 



lim X = a, 

r — >+oo 

where x = {f^+,R, Z) and a = ^5,0, , considering the cases 1 < 7 < 2 and 7=1 

simultaneously. 

We begin by defining the subset D in theorem B.l by 

T,l+R'^ + Z^ < 1. 

We now verify the hypotheses Hi and H2. Firstly, let w : [ro,+oo) ^ D be any 
C°[to,oo) function. Since limT-^+oo M(t) = it follows immediately that 



lim g{w{T),T)= lim M{t) (Bj, RB^, ZB^) 



= 0, 



showing that Hi is satisfied. Secondly, H2 is satisfied since the variables R and 
Z are bounded for all r ^ tq with tq sufficiently large. Therefore, since 

lim x{t) = a 

T— >+CX) 

for all initial conditions x(to) in D (see (B.ll)), theorem B.l implies that 

lim x{t) = a (B.13) 

T— > + 00 

for all initial conditions x{tq) in D. 

Finally, we need to show that any initial condition x{to) = (S_|_, i?, 
M(to), iP{to) for the DE (2.17), subject to O > and (2.21), determines an initial 
condition x{to) in D for the DE (3.12), so that (B.13) is satisfied. 

Indeed, since lim^-^+oo V' = +00, we can without loss of generality restrict the 
initial condition V'(to) to be a multiple of tt. This requirement can be achieved by 
simply following the solution determined by the original initial condition until this 
condition is satisfied. It follows from this condition, in conjunction with (3.8) and the 
restriction O > applied to (2.19), that 

(f:l + R^ + z^)\ ^ (j:1 + r^ + z^)\ < i, 

\ T / IX — To V T / It — To 

SO that x{to) G D. 



Appendix B.3. The limit of M/R at late times 
In analogy to (3.8), we define a variable M by 

M = M (1 + iiJ^M sin 2t/>) . (B.14) 
It follows from (2.17) that the evolution equation for M is of the form 

M' = -{Q + 2f:+ + MB]^f)M, (B.15) 

where Bj^j is a bounded function for r sufficiently large. By using (3.9) we obtain 

={-l + h{x,M,i,))(^), (B.16) 



RJ ' ^ ' ' \R 
where 

h{x, M, ^) = I - 2(3 - 3S+ + MB* 
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and -B* is a bounded function for r sufficiently large. It follows from (3.6), (3.8), (3.10) 
and theorem 3.1 that ]mi^^j^^h{x, M = 0. Consequently, (B.16) implies that 

^=0(e(-W)r), if 1<7<2, (B.17) 

as r — > +00 for any 5 > Q. Therefore, Vivcit ^+00 M / R = 0, which implies that 
limT-^+oo M/R = 0, on account of (B.14) and (3.8). 



Appendix C. Derivation of the cisymptotic expansions (4.1) and (4.2) 

In this appendix we give details of the derivation of the asymptotic expansions (4.1) 
and (4.2) for S+, R, Z, M and as r ^ +00. As was shown in appendix B.2, all 
solutions x(t) of the DE (3.14) in the set S (defined by (B.7)) are future asymptotic 

to the equilibrium point (E+, R,Z) = ^5, 0, . Since this equilibrium point is non- 
hyperbolic, centre manifold theory is required in order to compute the decay rates of 
x{t). In addition, the analysis has to be split into two cases, 1 < 7 ^ 2 and 7 = 1, as 
described in appendices C.l and C.2. We refer the reader to Carr (1981, pp 1-13) for 
an introductory discussion of centre manifold theory. 

The next step is to show that the asymptotic expansions for x{t) arc compatible 
with the non-autonomous DE (3.12). In order to show this wc need an asymptotic 
decay rate for M, which we obtain as follows. Equation (B.17) and the boundedness 
of R imply that M = 0(e(~^/2+'')^), as r +00, for any 6 > Q. Using (B.14) we can 
conclude that 

M = (e^-^/^+^)-^ , (C.l) 

as T — > +00, for any 5 > 0. We then write the non-autonomous DE (3.12) in integral 
form 

x{t) = x{to) + [ [f{x{s)) + g{x{s), s)] ds, (C.2) 

•/to 

and can verify that the expansions for x(r) are compatible with (C.2) using (3.13) 
and (C.l). 



Appendix C.l. 1 < 7 < 2 

In order to put the DE (3.14) in the canonical form required for application of centre 
manifold theory, we begin by making the translation 

{±+,R,Z)={y+^,x,z+^y (C.3) 

We can now apply theorem 3 in Carr (1981) to obtain an approximation to the local 
centre manifold through the point {x, y, z) = (0, 0, 0) given by 

y = -x^ - 2x^ - f + 0{x^), 

Upon applying theorem 1 in Carr (1981), we find that the governing DE for the fiow 
on the centre manifold is given by 

x' = -2x^ - 4x^ - 32x'^ + 0{x^). (C.5) 
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To obtain the asymptotic form of x{t) from (C.5), we follow a procedure analogous 
to the one outlined in section 3.1 of Carr (1981) who considers a related problem. It 
follows that 



(C.6) 



as T +OC', where C, is a constant that depends on the initial conditions (which 
can be set to zero upon appropriate translation of r). We have also included the 
second-order term in x{t) so that the reader can compute the decay rates in section 4 
to higher order if desired. The decay rates for R and Z then follow from (C.3), 
(C.4) and (C.6) in conjunction with theorem 2 in Carr (1981). The compatibility 
check based on (C.l) and (C.2), in conjunction with the transformation (3.8), then 
lead to the given decay rates for S+, R and Z in (4.1). 

To obtain the asymptotic form of M in (4.1), we substitute the known asymptotic 
expansions for x{t) into the evolution equation (B.15) for M, obtaining 

M'=[-1 + It-' (l + 0(i^))]M, 

on account of (C.l). We can now solve this DE and use (B.14) to obtain the decay 
rate for M . 

Finally, to obtain the asymptotic expansion for f2, we define a variable f2 according 

to 

O = f2 (1 - \R^M sin 2V') . (C.7) 
It follows from (2.17) that the evolution equation for O is of the form 

n' = [2Q - (37 - 2) + MB^]Q., 

where B^ is a bounded function for r sufficiently large. Furthermore, using (3.9) we 
obtain 

iL)' = |_3h-l) + MB.|iL. 

We can now solve this DE and then use the decay rates for R and Z and (C.7) to 
obtain the asymptotic expansion for Q as stated in (4.2). 

Appendix C.2. 7 = 1 

Unlike the case 1 < 7 ^ 2 where the centre manifold is one-dimensional, when 
7 = 1 the centre manifold of the equilibrium point (S_|_,^, Z) = ^^,0, is two- 
dimensional. Although the analysis is potentially more complicated, it can be reduced 
to a one-dimensional centre manifold problem by considering the DE (B.3) on the 
invariant surfaces Cl/{RZ) = 2k (see (B.12)), where fc > is a constant. Since, the 
analysis parallels that in appendix C.l we only provide a brief outline of the proof. 
To begin, we use this first integral along with (B.5) to solve for S+ in terms of R and 
Z in a neighbourhood around 11+ = 5 to obtain 

S+ = (1 - 2kRZ -R? - Z^fl"^. (C.8) 

Substituting (C.8) into the R and Z evolution equations in (B.3) and then making 
the transformation 

R = X, 

(C.9) 
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puts the system into canonical form. We find that the centre manifold through the 
point {x, y) = (0, 0) is approximated by 

y = i^' + + 10fc')a=^ + f (2 + + 0{x^), (C.IO) 

the corresponding DE on the centre manifold is 

x' = -Vskx"^ + - 2) - 2V3kx'^ + 0{x^), (C.ll) 
and the asymptotic expansion for x{t) is given by 



x{t) = — 



1 + ^ ((1 - QCl) Inr + a) + 0(t-2 In^ r) 



(C.12) 



as T — > +0O, where C* is a constant that depends on the initial conditions. We note 
that the constant Cr is related to the constant k through Cr = l/(\/3^)- 

The asymptotic decay rates for R and Z follow from (C.9), (C.IO) and (C.12), in 
conjunction with theorem 2 in Carr (1981). The expansion for is obtained from 
(C.8). Dropping the hats as in appendix C.l yields the decay rates for S+, R and Z 
as stated in equation (4.1). The asymptotic form of M is obtained in an analogous 
fashion to that outlined in appendix C.l. Finally, the asymptotic expansion for O as 
stated in (4.2) is obtained directly from (2.19). 

Appendix D. The Weyl curvature tensor 

In this appendix we give an expression for the Wcyl curvature parameter W in terms 
of the Hubble-normalized variables S+, i?, Z, M and i/i. Let Eap and Hap be the 
components of the electric and magnetic parts of the Weyl tensor relative to the group 
invariant frame with eg = u. Since Eap and Hap are diagonal and trace-free (see WE, 
chapter 6, appendix), they each have two independent components. In analogy with 
(2.3) we define 

^+ = 5(^22 + ^33), = -5:75(^22 - ^33), 

1 1 (D-1) 

H+ = 2(^22 + W33), H- = ^;^(W22 — W33), 

where £ai3 and Hap are the dimensionless counterparts of Eap and Hap, defined by 

It follows from (3.2), (D.l) and (D.2) that 

W'^ = £l + £l+nl + n'i. (D.3) 

The desired expressions for £± and H± can be obtained directly from equations (6.36) 
and (6.37) in WE, using (2.16) in the present paper: 

£+ = 5]+(l + E+) + ii?2(l - 3cos2V') -Z'^- 'iM^Z'^, 

H+ = -^R^sm2il) + ^MY:+Z^, 

£- = ^ [sinV + 5M(l-2E+)cosV + iM2z2sinV] , ^^'^^ 
1R 

n- = — [-cos^- |MS+sin^- fM^Z^cosV] • 
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